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Gas in Variational Perturbation Theory Through Six Loops 
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We compute the shift of the transition temperature for a homogenous weakly interacting Bose gas 
in leading order in the scattering length a for given particle density n. Using variational perturbation 
theory through six loops in a classical three-dimensional scalar field theory, we obtain AT C /T C — 
1.25 ± 0.13 an 1 ' 3 , in agreement with recent Monte-Carlo results. 
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A dilute homogenous Bose gas with particle density n, 
where the scattering length a is small compared to the 
interparticle spacing ~ n" 1 / 3 is a fine example of how, 
under the right conditions, collective effects can generate 
strongly coupled modes from microscopic degrees of free- 
dom exhibiting non-zero, but arbitrarily weak interac- 
tions. These conditions are met when the temperature is 
close to the transition temperature for Bose-Einstein con- 
densation (BEC). Naive perturbation theory (PT) then 
breaks down for physical quantities sensitive to the col- 
lective long-wavelength modes. One of the basic ques- 
tions is about the nature and size of the shift of the BEC 
transition temperature due to the interactions. Although 
attempts at the problem have a long history [l[ , only 
the recent advent of experimental realizations of BEC in 
dilute gases have prompted considerable work to finally 
solve the problem both qualitatively and quantitati vely 

iaaiaiiiiaiiiiBiiiiiiiiiiaeiiiiiiiif. 

Here we treat the case of a homogenous gas as opposed 
to, e.g., the case of a harmonic trap [2l[ . 

The appropriate formal framework for the treatment 
of the many-particle Schrodinger equation describing a 
gas of identical spin-0 bosons is non-relativistic (3 + 1)- 
dimensional field theory. Being interested only in equi- 
librium quantities, we switch to imaginary time r = it 
and consider the Euclidean action 
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where m is the mass of the bosons and \x the chemical 
potential and we work in units where ks = h = 1. For the 
description of the physics involving the long-wavelength 
modes, the modes with non-zero Matsubara frequencies 
may be integrated out from the theory, leaving a three- 
dimensional classical field theory. This is, in general, a 
complicated procedure, but, as a first approximation, we 
may simply ignore the dependence of ip on r and obtain 
the three-dimensional effective Euclidean action 
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with the identifications ip = V mT(^ + ifo), ?"barc = 
— 2m/i, u — A8namT. For details on the matching be- 
yond this approximation and why, for our purposes, it 
is indeed justified to use J2J), see Ref. It is conve- 
nient to generalize the model J5J) to an O(N) field theory, 



where <f> = (</>i, . . . , 4>n), 4> 2 = 4> a 4>a- This allows to make 
contact with the exactly known large- TV result for AT e 
0. 

It has been argued convincingly 0, 0] that the 
shift of the transition temperature away from the non- 
interacting value 
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and that the constant c\ is independent of the non-zero 
Matsubara modes and can be computed within the above 
three-dimensional field theory (for the exact value of 
c' 2 and an approximate evaluation of c'2, see Ref. |15)). 

Apart from the fact that the leading dependence of 
AT C on a has only been resolved recently , attempts at 
the determination of c\ have produced results from —0.93 
01 to 4.7 0] as has been summarized for instance in Ref. 
(13| . see also Fig. [5] below. The various approaches have 
different weaknesses and we comment on some of them 
below. It appears that the most reliable results so far are 
obtained by Monte-Carlo simulations (MC) [l2l[T3| . 

Improving perturbation theory by a variational princi- 
ple is an old idea |22|. Earlier att emp ts using resummed 
PT to determine AT C HH E3, EM have employed the 
linear 5 expansion (LDE). These were criticized in Refs. 
pfjl |23| . where Kleinert's field theoretic variational per- 
turbation theory (VPT, see [2(| I27L l28ll and Chapters 5 
and 19 of the textbooks |2|j and |30(. respectively) was 
used through five loops. We extend this calculation to 
six loops with a different treatment of the one-loop term, 
which is absent here. Our first non-zero coefficient arises 
only at three loops. In consequence, we obtain a 30% 
larger result than Ref. [2(j ■ We comment on the issue at 
the end of the paragraph following Eq. © below. 

After the interaction is turned on, we can either tune 
the temperature by AT C or the particle density by An 
to bring the system back to the transition, keeping the 
respective other quantity fixed. Eq. © relates the two 
shifts in leading order [7|, 
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In other words, we may determine c\ by computing An 
at the critical point. An is related to the interaction- 
induced shift of (0 2 ) in the three-dimensional theory by 



An = A(^) = mTA(<p 2 ). 
Combining ©, JU, © and ©, ci is given by 
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and the remaining task is to compute the critical limit of 
A((f) 2 }/Nu in the three-dimensional theory. We have 
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[G(p) - G (p)}, 
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where G and Go are the interacting and the free propaga- 
tor, respectively. The theory J5J) is super-renormalizable 
and only fbare has to be tuned to cancel UV divergences 
(a further divergence in the free energy is of no interest 
to us here). A convenient renormalization scheme de- 
fines the renormalized quantity r — /"bare — ^(0), where 
S(p) is the self-energy. The full propagator reads then 
G(p) = l/{p 2 +r-[E(p)-S(0)]} and r = corresponds to 
the critical point, where BEC takes place. Consequently, 
Go(p) = l/{p 2 + r) 1191 . A nother pop ular scheme is to 
use G (p) = 1/p 2 jnTHl E ES HI. This generates 
an unnatural one-loop contribution to AT C even in the 
absence of interactions. Although this contribution van- 
ishes as r — > 0, it strongly influences rcsummation and is 
responsible for the small value of ci in Ref . (2£j • 

Perturbation theory organizes A(<fi 2 )/Nu in a loop ex- 
pansion 
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At one loop, there is no interaction and therefore the 
difference in JSJ vanishes. There is also no two-loop 
contribution since our renormalization scheme subtracts 
all graphs containing momentum-independent self-energy 
insertions and there is no two-loop graph without such an 
insertion. This reflects the fact that the one-loop correc- 
tion of the propagator merely renormalizes the chemical 
potential, which does not lead to a shift in the critical 
temperature |7J, llg • The three-loop contribution is 
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where the dashed box means that the enclosed self-energy 
contribution has to be taken at zero external momentum. 
The three-loop case exhibits the only type of divergence 
appearing at any loop order of the calculation, the "sun- 
set" subdivergence which causes the two diagrams in 
to be UV divergent. One may either regulate each dia- 
gram and remove the regulator from the finite difference 
or, much simpler, perform the momentum integrations 
only after the subtraction to arrive at the given result. 



To save space, it is convenient to define an operator 
1Z that recursively subtracts the zero-momentum contri- 
butions of all self-energy insertions of a given diagram. 
E.g., the combination of diagrams in Eq. IjlOl) is then rep- 
resented by 1Z applied to the first diagram. The diagrams 
through six loops, their weights and 0(A) symmetry fac- 
tors as well as the numerical values of the corresponding 
integrals are given in Table H The latter were obtained 
by own computations and confirmed by extracting them 
from Ref. [2J|. Details will be presented elsewhere [25| . 
The values for the diagrams through five loops are consis- 
tent with those in Ref. 01 [ we have more precise results 
for diagrams 5-4 and 5-5, though; the five-loop results 
of have to be taken with care since at least two of 
the diagrams (5-2 and 5-3) are wrong there, while the 
corresponding coefficient as is off by 10%]. 

The resulting perturbative coefficients are a\ — a 2 = 
0, a 3 = -5.06047 x 10" 5 , a 4 = 2.99626 x 10~ 6 , a 5 = 
-1.93583 x 10~ 7 , a 6 = 1.31373 x 10~ 8 for N -> oo and 
ai = a 2 = 0, a 3 = -1.01209 x 10~ 4 , a 4 = 2.99626 x 10~ 5 , 
a 5 = -1.19872 x 10~ 5 , a 6 = 5.85519 x 10" 6 for N = 2. 

To extract the r — > limit of @, we need to apply 
a resummation algorithm. The alternating signs of the 
coefficients a; suggest to use a Borel-type method, but 
since we have no information on the distance of the first 
singularity of the Borel sum from the origin in the com- 
plex plane, we proceed differently and use VPT. This 
method has proven to yield accurate results for critical 
exponents [27l I28L Effll ] and amplitude ratios m the 
present context of critical phenomena. The crucial at- 
tractive feature of VPT for our purpose is that it allows 
to implement the correct Wegner exponent of approach to 
scaling to [— (3'(g*) in a renormalization group treatment] 
[22]. For quantities like A((j) 2 )/Nu that remain finite at 
the critical point, this approach has, in our renormaliza- 
tion scheme, an expansion 
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where lu' = 2lu/(2 — n) and r\ is the anomalous dimension 
of the critical propagator. 



For a truncated partial sum Yli=i a i 
Klcinert's VPT requires replacing 

{tu) l - 2 \l + t ( — ] I 
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(note that this is an identity for t = 1), reexpanding 
the resulting expression in t through t L ~ 2 , setting t = 1 
and then optimizing in u, where optimizing is done in 
accordance with the principle of minimal sensitivity [331 ] 
and in practice means finding appropriate stationary or 
turning points. That is, we replace 
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and optimize the resulting expression in u. For u r — > 00, 
we obtain the L-loop approximation of the r — ► limit 
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iV(iV+2) 
3 



-7.24858 x 10" 



iV(iV+2)(iV+8) 
27 



-2.04919 x 10" 6 
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2.32536 x 10" 





6-1 
ft 




6-2 

n 



6-3 

n 



6-4 
ft 



6-5 
ft 






1 iV(iV+2)(5iV+22) 
4 81 



-3.47299 x 10" 



1 iV(iV+2)(5iV+22) 
4 81 



-5.39050 x 10" 



iV(iV+2)(iV 3 +8iV 2 +24JV+48) 
243 



-2.04830 x 10 -9 



N(N+2Y(N+S) 
81 



1.83919 x 10" 



N(N+2Y (N+S) 
81 



8.98169 x 10" 



N(N+2Y(N+S) 
81 



1.40475 x 10" 



N(N+2)' ! (N+8) 
81 



1.23133 x 10" 
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iV(iV+2)(3iV z +22JV+56) 
243 



-1.10895 x 10" 9 



N(N+2)(3N' ! +22N+56) 
243 



-1.65022 x 10" 



iV(iV+2)(3iV z +22JV+56) 
243 



.05795 x 10" 



N(N+2)(N^ +20iV+(i()) 
243 



-6.27851 x 10" 



N(N+2)(N^ +20iV+(i()) 
243 



-1.06164 x 10" 



N(N+2)(N' ! +20N+6{)) 
243 



-1.32467 x 10" 



i_ N(N+2)(5N+22) 
6 81 



-4.25366 x 10" 



TABLE I: Diagrams from three through six loops, their weights wl-u, O(N) symmetry factors gL-n and numerical results II- 
for the corresponding integrals for r = 1. The contribution of each diagram to (4> 2 ) is (— u) L ~ 
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Results are only available starting at four loops, since 
two non-zero perturbative coefficients are necessary for 
VPT to work. 

We first treat the N — > oo limit, where rj = and 
lu = lu' = 1 (see, e.g., Ref. 0). For large N, the only 
surviving diagrams are 3-1, 4-1, 5-1, 6-1 and their gen- 
eralization to higher orders, which are all easy to com- 
pute. The corresponding resummed results start out 



as /< 4) = -t 

p(6) _ 



.54677 x 1(T 4 , = -9.14163 x 10" 4 , 

F (40) 



/^ D; = -9.55894 x 1(T 4 . We have plotted /^ 4) , . . . , 
in Fig.Q] They are seen to converge, albeit slowly, to the 
exact result A(0 2 )/7Vm = -l/96vr 2 = -1.05543 x 1CT 3 
@. We will return to the issue of slow convergence for 
the iV — > oo limit below. 



arbitrarily setting u>' = 2. The corresponding values start 
out as /g 4 ^ = 



f(6) 
JO 



-4.80756 x 10- 4 , /($ 5) = -5.53705 x 10 



-6.11288 x 10~ 4 . We have included 



(4) 



• ./ 



(40) 



in Fig. n to demonstrate the at best very slow conver- 
gence. There are attempts to accelerate it 0, but it 
has been argued in Ref. (23 that LDE is inapplicable to 
field theory whenever there is an irrational exponent of 
approach to scaling. 

Now we turn to the case N = 2. We employ two ap- 
proaches when determining the r — > limit of A(</> 2 ) /Nu. 
One of them uses the known values u = 0.79 ± 0.01, 
77 = 0.037 ± 0.003 (see, e.g., Ref. [U) and thus w' = 
0.805 ± 0.011 for N = 2, and the other self-consistently 
determines us' at each order. With ui' = 0.805 (its un- 
certainty ±0.011 causes a 1% uncertainty of the fixed-a/ 
six-loop estimate for c\) we obtain /q 4 -* = —4.297 x 10~ 4 , 
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lation to c\ via 10 is plotted in Fig. [2 



.103x10 4 , whose trans- 
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FIG. 1: Resummed results for A{cj> 2 ) /Nu as a function of 
the number of loops L for TV — > 00. Upper dots: LDE. Lower 
dots: VPT. Solid line: exact result. 

The LDE used in Refs. [H 13 El E3 corresponds to 
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FIG. 2: ci as a function of the number of loops L for N = 2. 
Upper dots: self-consistent uj' . Lower dots: ui' = 0.805. The 
label VPT' indicates the inclusion of a non-zero one-loop term 
in Ref. |2fl ]. which is absent in the present treatment, labeled 
by VPT. Also included are most results from other sources. 
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The expected large-L behavior in VPT has the form 

/ +Acxp[-.B(L-3) 1 - cy ] IHHIS. However > even and 
odd loop orders in VPT are, in our case, obtained from 
extrema and turning points, respectively, which causes 
them to behave slightly differently. This prohibits a 
trustworthy extrapolation to L — > oo at the current loop 
order. Note that in the large- N limit, where u>' = 1, the 
expected large-L behavior turns into fo + A' /(L — 3) s 
and an analysis of the data shown in Fig.^ but extended 
through /^ 150) , yields A 1 w 3.4 x 1CT 4 , B' « 1.0. 

For determining u>' self-consistently [53, H3 , we apply 
VPT to the quantity dln(A((/> 2 ) /Nu) /d\au r and exploit 
the fact that it vanishes as u r — > oo and that its large-it,, 
expansion has the same w'-exponent as A(</> 2 ) /Nu. That 
is, at each loop order we tune to' such that the optimiza- 
tion procedure leads to a vanishing approximation for 
the large-u r limit of dln(A(0 2 ) /Nu) / d\nu r . This u>' is 
then used at the same loop order in ((Hfl . Our results are 
f^ 5) = -6.343 x 10- 4 and = -6.268 x 10" 4 , with cor- 
responding w'-values = 0.6212 and uj'^ = 0.6381. 



Via © they are converted to the estimates for c\ plotted 
in Fig. 

Taking the average and the difference of our two esti- 
mates for /q 6 ^ as the mean value and the total error bar 
of our final result, respectively, we arrive at A(0 2 ) /Nu — 
-(5.69±0.58) x 10~ 4 and via at c x = 1.25±0.13. This 
is in agreement with the recent MC data ci = 1.32 ±0.02 
[13 and ci = 1.29 ± 0.05 (see Refs. EH Q for ex- 
planations of the low MC value in Ref. see Ref. 
for a critique of the MC calculation of Ref. Q). Fig. [21 
shows a comparison with most available results. 

In summary, we provide a theoretically well-founded 
resummed perturbative result for c%, comparing favor- 
ably with recent Monte-Carlo data. 

The author is grateful to H. Kleinert for important 
suggestions, many discussions, and a careful reading of 
the manuscript; also to E. Braaten for many helpful com- 
ments on the manuscript and to B. Nickel and T. Binoth 
for useful communications. 
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